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Abstract 

Given uj > 1, let Z? w j be the graph with vertex set 1? in which two 
vertices are joined if they agree in one coordinate and differ by at most ui in 
the other. (Thus Zj^v is precisely Z 2 .) Let p c (u) be the critical probability 
for site percolation on Z 2 u y Extending recent results of Frieze, Kleinberg, 
Ravi and Debany, we show that lim^^oo u)p c {uj) — log(3/2). We also prove 
analogues of this result for the n-by-n grid and in higher dimensions, the 
latter involving interesting connections to Gilbert's continuum percolation 
model. To prove our results, we explore the component of the origin 
in a certain non-standard way, and show that this exploration is well 
approximated by a certain branching random walk. 



1 Introduction and results 

Frieze, Kleinberg, Ravi and Debany [9] proposed the following random graph 
as a model of an ad hoc network in an environment with (regular) obstructions. 
Given positive integers n and u, let G — [n]?^ be the graph with vertex set 
[n] x [n] in which two vertices {x\,yi) and (2:2,2/2) are joined if X\ — X2 and 
1 2/x — 2/2 1 < w or Di — Vi an d \x\ —xi\ < uj. Let V be a random subset of [n] x [n] 
obtained by selecting each point (x, y) with probability p, independently of the 
other points, and let G[V] denote the subgraph of G induced by the vertices in 
V. We write G n ,uj, P for G[V], and G UtP for the infinite random subgraph of 
defined in the same way but starting from Z 2 rather than [n] x [n] . 

The model G' n uj p studied in [S] was defined on the n by n torus rather than 
the grid, but G n>ul ,p seems more natural in terms of the original motivation. For 
the results here, it will make no difference which variant we consider. 

One interpretation of G n ,uj, P is as follows: sensors are dropped onto a random 
selection of crossroads in a regularly laid out city (or planted forest); two sensors 
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can communicate if they are within distance to and the line of sight between them 
is not blocked by a building (tree). The graph G n ^, P indicates which pairs of 
sensors can communicate directly, so we would like to know when G„ tUJ ,p has a 
giant component, and roughly how large it is. Alternatively, taking the point 
of view of percolation theory, we would like to study percolation in G^ iP , i.e., 
to know for which choices of the parameters G WiP has an infinite component. 

Kolmogorov's 0/1-law implies that for each fixed u> there is a 'critical prob- 
ability' p c {u) such that if p > p c (oj) then G WiP contains an infinite component 
with probability 1, while if p < p c (uj) then with probability 1 all components of 
G w , p are finite. Furthermore, by the 'uniqueness theorem' of Aizenman, Kesten 
and Newman [T] (see also p. 121]), when an infinite cluster exists it is unique 
with probability 1. Among other results, Frieze, Kleinberg, Ravi and Debany [9] 
proved (essentially) that 

l/(4e) < liminf uop c (uj) < limsupu>p c (w) < 1, 

and posed the natural question of determining the value of lmiujp c (uj), assum- 
ing it exists. (In fact, they worked with finite random graphs, as we shall in 
Theorem [5] below.) Here we shall answer this question, proving that 

lim lo Pc {lu) = log(3/2) = 0.4054 . . . , (1) 

u — >oc 

as part of the more detailed Theorem Q] below. 

In homogeneous cases such as this, where the expected degrees of all vertices 
are equal, the critical probability is perhaps more naturally described in terms 
of the 'critical expected degree'. Since the expected degree of a vertex of G U)P 
is 4pw, jl]) states that the critical expected deg 166 of Guj^p tends to 41og(3/2) = 
1.6218 ... as w -* oo. 

Note that it is too much to hope to find p c {io) exactly for a given lo: with 
lo fixed, Pc(ijj) is the critical probability for site percolation on a certain lattice 
^(o)) (non-planar for ui > 2); the exact values of such critical probabilities are 
known only in a few very special cases. Indeed, p c (l) is the critical probability 
for site percolation on Z 2 , and even this is not known exactly. 

In spirit, Theorem[T]and its proof are very close to the arguments given in [BJ 
for Penrose's spread-out percolation model 13J : when lo — > oo withpcj bounded, 
we can approximate the neighbourhood exploration process in G UjP by a certain 
branching process, and describe the distribution of the small components in 
this graph in terms of the branching process. We then show by a rather ad 
hoc argument that almost all vertices not in 'small' components are in a single 
infinite component. An analogous assertion holds for G n ,uj, P and its unique 
giant component. In contrast to [BJ, where we could appeal to the general 
sparse inhomogeneous random graph model of [5], here even the local coupling 
to the branching process has to be done by hand. We shall discuss this later. 

To state our results precisely, we first define the branching process alluded 
to above. Let /i > be a real parameter, and let T = T M denote (a random 
variable with) the geometric distribution with parameter 1 — e~''. Thus, 

P(r = k) = (1 - e-^fe-* 1 
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for k = 0, 1, . . ., and E(T) = e' 1 — 1. Let denote the distribution of the sum 
of r independent copies of T M . 

Let Xfj, be the branching process (Xq, Xi, . . .) in which Xq consists of a 
single particle xq, each particle x € X t has children independently of the other 
particles and of the history, the number of children of xq has the distribution 

, and the number of children of a particle x £ X t , t > 1, has the distribution 
rjj . Let <fi([i) denote the survival probability of the branching process 3E„, i.e., 

00) = ^{\Xt\ > for all t). 

Note that 

E(rf) = 2E(F M ) - 2(e^ - 1), 

so the branching process X^ is supercritical if and only if fi > log(3/2), i.e., 
4>{p) > if and only if /i > log(3/2). 

In the result below, Co denotes the component of the origin in the random 
graph G WiP ; if the origin is not a vertex of this graph, then we set Co = 0. 
We write P w , p for the probability measure associated to G w , p , and 9(uj,p) — 
^uj,p(\Co\ = oo) for the probability that Co is infinite. 

Theorem 1. Let < A < log(3/2) be constant. If uj is large enough then 
all components of G^.x/u finite with probability 1. Furthermore, there is a 
constant a\ > such that, for all large enough u, 

^.x/U\Co\>k)<-e- a ^ k -^ 

holds for every k > 1. 

Let A > log(3/2) be constant. If uj is large enough, then G u ^\/ U contains a 
unique infinite component with probability 1. Furthermore, 

6(u,\/uj)~<f>(X)\/u; (2) 

as uj —> oo. 

Note that Theorem [T] certainly implies (fi} , which is equivalent (as noted 
earlier) to the statement that the critical expected degree tends to 41og(3/2) 
as uj — > oo. The proof of Theorem [T] will form the bulk of Section [3J In fact, 
the upper bound on component sizes (giving the lower bound on the critical 
expected degree), is fairly easy. The lower bound on component sizes (giving 
the upper bound on the critical expected degree) will be derived in a way that 
will enable us to prove a finite analogue of Theorem [TJ stated as Theorem [2] 
below. 

For the upper bound on component sizes we shall prove a much more detailed 
statement than stated in Theorem [T] namely that 

P UJ ,(|Co|>*)<pP(|3£ M |>fc) (3) 

for every w>l, 0<J»<1 and k > 1, where 

e-" = (l- P r, (4) 
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and |X M | denotes the total number of particles in all generations of 3£ M . For 
/i < log(3/2), the branching process is strictly subcritical, and standard 
results imply that there is some 6 M > such that P(|X M | > k) < e~ fcf '( fe_1 ' holds 
for every k. If A < log(3/2) is constant and we set p = X/uj, then as u> — > oo 
we have [i ~ pu> — A. In particular, /i is bounded above by some // < log(3/2) 
when u is large enough, so the first statement in Theorem [1] follows. 

Letting k tend to infinity in ([3]) (or simply setting k = oo, which makes 
perfect sense), we see that 6{uj,p) < ptf>(ii). Standard results on branching 
processes imply that (j>(-) is continuous. Thus ([3]) implies the upper bound 
implicit in @, i.e., that 

6(w, X/uj) < (1 + o(l))0(A)A/w 

as — > oo with A > log(3/2) constant. 

Note that for the exact comparison with a branching process above, it is 
natural to work with the parameter fi = — log((l — p) u ). However, in the lower 
bounds on component sizes, we can obtain only approximate results, and it is 
more natural to use the asymptotically equivalent parameter A = pu), which we 
have used throughout the statement of Theorem [TJ We shall switch between 
these two parameters as and when convenient. 

Although our main focus is the infinite random graph G WjP , our other aim 
is to prove a result corresponding to Theorem [T] for the finite graphs G n)W)P 
and G' u . Note that there is a natural coupling in which G n>UJtV , defined on 
the grid, is a subgraph of G' nup , defined on the torus. Let C\(G) denote the 
number of vertices in a largest component of a graph G. 

Theorem 2. Let < A < log(3/2) be constant. There are constants A\ and 
uj\ such that 

Ci(G n<UtX/ui ) < Ct{G' niUtX/ J < A x logn 

holds with probability 1 — o(l) as n ^ oo, whenever w = U)(n) satisfies u>\ < 
lu < n. 

Let X > log(3/2) be constant. If uj — u>(n) is such that lo and n/u) — > oo, 
then 

Ci(G„,«,A/a,), ^(G'^J = (0(A) + o p {l))nX/u 

as n — > oo. 

We shall also state results for the case uj — > oo with n = Q(u>); for details 
see Subsection 12.31 

2 Proofs 

2.1 Upper bounds on component sizes 

We start with some simple upper bounds on the component sizes. Let p and 
uj be arbitrary. To prove ([3]), it suffices to couple G UtP and X M so that if we 
explore the component Co of the origin in G UtP in a suitable manner, then the 
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number of vertices reached at each stage is at most the number of particles of 
X M in the corresponding generation, where /i = — wlog(l — p) is given by (|4]). 

In each step of our overall exploration process (except for the first) we either 
'explore vertically', first up and then down, or 'explore horizontally', first to the 
left and then to the right. To 'explore upwards' from a vertex v = (x, y) of G UtP , 
we test the points (x,y + 1), (x,y + 2), . . . , one by one. Each test 'succeeds' 
if the relevant point is a vertex of G u .p that was not previously reached (in an 
earlier stage of the overall exploration). We stop the upwards exploration as 
soon as ui consecutive tests fail. 

If Ui, . . . , v r denote the vertices of G WjP reached during this upwards explo- 
ration, corresponding to the successful tests, then vv\, V1V2, ■ ■ ■ , v r -\v r are all 
edges of G W;P . Also, there is typically no edge of G WiP from v r — [x,y') to a 
point (x, y") with y" > y'; this is not always true as a test may fail due to the 
presence of a vertex we reached earlier. If we have not previously tested the rel- 
evant points, so each is present in G WiP with probability p, independently of the 
others, then the probability that the next u> tests fail is exactly (1 — p) u = e _At , 
so r has the geometric distribution T^. In general, we may have previously 
tested some of the points; a repeated test of a point always fails by definition, 
so, conditional on previous steps of the exploration, the distribution of r is 
stochastically dominated by T^. 

To 'explore vertically' from a vertex v, we explore upwards from v, and then 
explore downwards from v. Similarly, to explore horizontally, we first explore 
to the left, and then to the right. 

Let us condition on vq — (0, 0) being a vertex of G^p, an event of probability 
p. Set Yq — {vq}- Explore vertically and horizontally from vo, and let Y\ be the 
set of vertices of G u . p reached in these explorations. From the remarks above, 
|Yi| has exactly the distribution Tjf ' . Suppose that we have defined Y u t > 0. 
For every vertex v G Y t that was reached during a horizontal exploration, we 
explore vertically from v. Similarly, we explore horizontally from each v G Y t 
reached during a vertical exploration. Let Yt+i be the set of new vertices of 
G w , p reached. 

It is not hard to check that the exploration just defined does indeed uncover 
the whole component Go- Indeed, let C = Ut>o^- Then we certainly have 
vo G Gq C Go- Since Go is a connected subgraph of G w ,p, if G ^ Go then there 
are vertices v G C' and w G Go \ C' that are adjacent in G^.p- Choosing v and 
w at minimal (Euclidean) distance, there are no vertices of G WiP between v and 
w. We may suppose without loss of generality that v = {x,y) and w = {x',y) 
with < x' — x < lu. When exploring, we reach v at some stage; if we do so 
during a vertical exploration, then exploring from v to the right we reach w at 
the first step. On the other hand, if we reach v while exploring to the right or 
to the left, the same exploration would also reach w, either just after or just 
before reaching v. This shows that w G G , contradicting our assumptions. 

From the remarks above, conditional on previous explorations, the number 
of vertices reached in a horizontal or vertical exploration is stochastically dom- 
mated by the distribution. Thus the whole sequence |1q|, |Yi|, |F 2 |, . . . is 
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stochastically dominated by \X \, \Xi\, \X 2 \, ■ ■ ■ , where X ,Xi, . . . are the gen- 
erations of the branching process X M . Since Co = U*>o we thus have 

P(|C | >k\v e V(G U , P )) < P(|XJ > k) 

for every k > 1, which gives ([3]). 

As noted earlier, the first part of Theorem [T] follows, along with the upper 
bound implicit in @. 

2.2 The lower bound on component sizes 

Throughout this section we shall set p = X/oj, where A > log(3/2) is constant 
and we let lu — > oo. All asymptotic notation refers to this limit. 

Usually, in showing that the neighbourhood exploration process in some 
random graph is well approximated by a branching process, one first shows that 
a given vertex is unlikely to be in a short cycle. Here the latter statement 
is false: given that v = (x,y) is a vertex of G^p, it has probability 6(1) of 
being in a triangle of the form (x,y), [x,y + a), (x, y + b), < a < b < lu. 
However, these triangles (and similar cycles) have already been accounted for 
in the construction of the branching process. If we delete all edges of G^ iP 
corresponding to line segments whose interiors contain vertices of G UtP , then it 
is easy to check that the density of short cycles in the resulting graph is very 
low. We shall show something more or less equivalent to this statement. 

Let k be fixed. As before, let us condition on the event that i>o = (0, 0) is 
a vertex of G UtP . We claim that we can couple the first k steps Yq, Y\, . . . ,Yk 
of the exploration defined above with the first k generations of X\ so that they 
agree (in the sense that \Y t \ = \X t \ for < t < k) with probability 1 — o(l). 
(Here the asymptotic notation refers to the limit lu — > oo with k and A = puj 
fixed.) We have already shown that the exploration is stochastically dominated 
by X M , where /i = — o;log(l — p). Since [i ~ A, the processes X M and Xa may be 
coupled so that their first k generations agree whp, so it suffices to 'bound (Y t ) 
from below', i.e., to couple a subset of this exploration process with Xa so as to 
agree whp for k generations. 

In doing so, we must be careful: suppose that as part of our exploration we 
find vertices of Co located as in Figure [1] Then we are in trouble - since w and 
w' lie on the same horizontal line, the horizontal explorations from these points 
will necessarily interfere. In particular, the expected number of new vertices 

(2) 

reached from w and w together will be significantly less than 2E(r^ ). Thus, 
if we reach w first, we do not allow ourselves to test the point w' . 

In general, we consider the lines of Z 2 , i.e., the sets of the form {x} x Z or 
Z x {y}. Initially, we mark all lines as 'unclaimed', except the two lines through 
v , which we mark as 'claimed by vq\ Whenever we find a vertex v £ Co during 
a horizontal exploration, this vertex immediately claims the vertical line it lies 
on. Similarly, when we find a vertex v during a vertical exploration, v claims the 
horizontal line it lies on. We modify our exploration as follows: when exploring 
horizontally or vertically from a vertex w, we omit testing any point w on a 
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Figure 1: Part of our exploration, including a possible clash (between w and 
w') that we must avoid. The solid lines indicate regions that have already been 
tested for vertices of G^ iP ; these lines extend a distance u> from the last vertex 
found. 



line already claimed by a vertex other than v. This ensures that the same line 
cannot be claimed twice: when we find a vertex w G C by exploring from v, 
exactly one line through w has been claimed, namely that joining v to w; the 
vertex w then claims the other unclaimed line. 

In this modified exploration, a point (x, y) can only be tested once, so each 
tested point is present in G WjP with (conditional) probability p. Since the ex- 
pected number of vertices reached within k steps of the exploration is 0(1), whp 
there are at most log to such vertices, say, and hence at most 2 log w claimed lines. 
In each horizontal or vertical exploration, we are exploring along a claimed line, 
and we omit testing at most one point on each other claimed line. Thus we 
omit testing O(logw) = o(l/p) points. Since we would have expected only o(l) 
successes among these tests, omitting them makes essentially no difference: the 
number r of new vertices found may be coupled with a T y x ' distribution so as 
to agree with probability 1 — o(l). As we expect to reach O(l) vertices during 
the first k steps of our exploration, the sum of these o(l) error probabilities is 
still o(l), so the claim follows. 

The truth of the statement ' | C 1 > fc' is certainly determined by the numbers 
\Y \, . . . , |Vfc_i|. Thus, for every 1 < k < oo we have 

F.,a/4I C o| >k\v e V(G U>X/U )) - P(|X A | > k) 

as lo — > oo, i.e., 

^,X/U\C \ >fc)~P(|*A| >k)\/uj. 

It follows that there is some K = K(u) —> oo such that 

P w ,A/«-(|Co| >K)~F(\X X \ >K)X/w. 
Since the branching process 3£\ does not depend on ui, the right-hand side above 
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is asymptotically <p(X)X/u). To complete the proof of |(3J| it remains to show that 
K,x/ U {K < \C \ < oo) = o(A/w). 

In other words, roughly speaking, we must show that almost all vertices of 
G^j, x/ui m 'large' components are in infinite components. As noted in the in- 
troduction, the probability that there is an infinite component is either or 1, 
and when there is an infinite component it is unique with probability 1, so the 
remaining assertion of Theorem [T] follows. 

So far we have coupled the neighbourhood exploration process in G UiP with a 
branching process X\. If we keep track of the locations of the points in Yq, Y\, . . . 
as well as their number, then the appropriate limit object is a branching random 
walk. 

Indeed, let us turn Xa into a branching random walk Xa on R 2 (or, more 
formally, on R 2 x {v, h}) as follows. Each generation X t , t > 1, will consist of a 
finite set of points of R 2 , each labelled with either V (for reached by a vertical 
step) or 'h' (reached by a horizontal step). A particle of type 'v' at a point (x, y) 
has children in the next generation according to the following rule: generate a 
Poisson process V of intensity AonRx {y}. Starting from (x,y) and working 
to the right, include as children of (x, y) all points of V until the first time that 
we come to a gap of length greater than 1. Do the same to the left. All children 
have type 'h'. The rule for children of particles of type 'h' is similar, using the 
line {x} x R instead. Sometimes we start with two particles at the origin, one 
of type 'h' and one of type 'v'; we write X| for this branching random walk. 
At other times, we consider the same branching rule but start with a single 
particle, writing X\. When the starting rule is clear or unimportant, we write 
simply X\. Since the number of offspring of a particle in Xa is independent of 
the position of this particle, we may view X| as our original branching process 
Xa by simply ignoring the positions of the particles. 

It is easy to check that the rules above correspond to a certain limit of the 
horizontal and vertical explorations we considered earlier, where we take ui —* oo 
and rescale by dividing the coordinates of our lattice points by oj. In particular, 
the coupling argument above shows that for any fixed k, we may couple (Y t )t<k 
with the first k generations (X t )t<k of the branching random walk X 2 so that 
with probability 1 — o(l) there is exactly one particle (x',y') £ X t for each 
(x, y) £ Y t , 1 < t < k, and \x/ui — x'|, \yjio — y'\ < k/uj. 

We shall only study the branching random walk Xa in a trivial way, calcu- 
lating certain expectations; moreover, we shall avoid all detailed calculations, 
needing only 'soft' arguments. 

We start with a simple lemma showing that the supercritical branching ran- 
dom walk X M , fi > log(3/2), remains supercritical even when restricted to a large 
enough square. By the restriction of X M to a region R we mean the branching 
random walk obtained from X M by deleting all particles that lie outside R, along 
with all their descendants. We write X /i (C) for the restriction of X M to [— C, C] 2 . 
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Lemma 3. Let [i > log(3/2) and e > be fixed. There are constants A, B and 
T\ (depending on [i) with the following properties. 

(i) For any point (xq, jjq) of [—A, A] 2 , if we start the branching random walk 
X M with a single point ( of type 'h ', say) at {xq, yo) and restrict to [-B, B] 2 , then 
the expected number of points of generation T\ inside [—A, A] 2 is at least 2. 

(ii) We have 

P(X^(-B) survives forever) > (p([i) — 3e, (5) 
where </>(//) is the survival probability of the branching process X M . 

Proof. Starting the unrestricted walk X* with a single particle at the origin, the 

expected size of generation t is E(r| 1 2 ^)' = (2e M — 2)*. Since fj, > log(3/2), there 
is some T\ such that this expectation is at least 10. Fix such a T\. 

Let X* (^4) denote the restriction of X^ to [—A, A] 2 . As A — > oo, the first T\ 
generations of X* (^4) converge in distribution to the first T\ generations of X* . 

Hence there is some Aq such that the expected size of generation T\ of X^(A) 
is at least 8 whenever A > Aq. 

By symmetry, the distributions of the number of points of generation T\ 
of X^(j4o) m each of the four squares [0, ±^4o] x [0, ±^4o] are identical. Let V 
denote this distribution, noting that E(X>) > 2. 

Let A > Ao be a constant depending on n and e, to be chosen later. Set 
B = 2A, and let (xo,yo) be any point of [— A, A] 2 . Without loss of generality, 
we may assume that Xo,yo < 0, so [xq,Xq + A] x [yo,yo + A] is contained in 
[— A, A] 2 ; see Figure [5] Let X^ denote the branching random walk X* started 
at (xo,Vo) an d restricted to [— B,B] 2 . Since the square of side 2A centered at 
(^0)2/o) is contained in [— B,B] 2 , the process X^ stochastically dominates the 
translated process X*(.A) + (xo,yo). In particular, the number of points of X^ 
in [—A, A] 2 , which is at least the number of points in [xo, xq + A] x [y , y + A], 
stochastically dominates V. Since E(D) > 2, this proves (i). 

We have shown that, for any A > Aq, if we start our branching walk re- 
stricted to [— 2A,2A] 2 at any point of [— A, A] 2 , and look at the number of 
points of [—A, A] 2 reached at times 0, T\, 2Ti, 3Ti,. . . , then the numbers we see 
stochastically dominate a certain supercritical Galton-Watson branching pro- 
cess (with offspring distribution T>). This process survives forever with some 
probability po > depending on fj, only. 

Recall that ignoring the positions of the particles turns the (unrestricted) 
branching walk X^ into our original branching process X M . Except that the off- 
spring distribution is different for the first generation, X^ is a standard Galton- 
Watson process, so standard results (for example, Theorem 2 in Athreya and 
Ney [2]) imply that with probability 1, either X M dies out (an event of proba- 
bility 1 — <p(n)), or \Xt\ — > oo as t — > oo. It follows that there is a time T such 
that with probability at least cf)([i)—£ generation T contains at least log(l/e)/po 
particles. Hence, choosing A > Aq large enough, the walk X^(^4) has probability 
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2B 



Figure 2: The dotted lines show the square S of side 2 A centered at (xo,yo), 
which is contained within the outer square of side 2B — 4A centered at the 
origin. The shaded quadrant of S is contained in [—A, A] 2 . 



at least ^(m) ~~ 2e of generating at least log(l/e)/po particles in generation T; 
these particles lie in [—A, A] 2 by definition. 

Let iV denote the number of particles of X^B) in generation T in [—A, A} 2 . 
Since X^B) may be regarded as a superset of X^A), we have N > log(l/e)/po 
with probability at least </>(/i) — 2e. Each particle of X'i(B) in [—A, A] 2 survives 
forever with probability at least po, independently of the others, so 

P(X 2 (B) dies out) < P(N < log(l/e)/ Po ) + (1 - p ) los(1/e)/po 
< l-(0(f*)-2e)+e=l-(0(fi)-3e). 

In other words, X^B) survives with probability at least </>(//) — 3e, proving 
©. □ 

In the next lemma we consider the restricted walk X^B) started at a particle 
v of type 'h' just below the bottom edge of [— B,B] 2 ; making sense of this 
requires a slight modification of our definition of restriction. Recalling that 
we label particles by the kind of exploration that reached them, the first step 
from v is vertical, so the first generation X\ of the unrestricted walk started at v 
consists of points on the same vertical line as v. To define the first generation X[ 
of our restriction X l _ l (B), first delete all points of X\ outside [— B,B] 2 . Then, 
if some point in the remaining set Xi D [— B,B] 2 is within distance 1 of v, 
let X[ = Xi n [-B,B] 2 ; otherwise, delete all remaining points too and let 
X[ = 0. All points in X[ are inside [— B,B] 2 , and from here we continue the 
restricted walk as usual. The reason for the somewhat fussy definition of the 
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first generation is that when we return to the graph G u>p , we will wish to find 
certain paths starting at a vertex v just outside [—u)B, toB] 2 , with all remaining 
vertices inside [— loB, wB] 2 ; the first edge of such a path must join v to a point 
in [— loB, u>B] 2 at distance at most u> from v. 

T * Si 
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Figure 3: The shaded regions are the sets Sq and Si considered in Lemma [4] 
and Corollary [5j the solid square is the boundary of [— B,B] 2 . The paths are 
for T2 — 16. Note that there may be additional points of the branching process 
in the interiors of the line segments of these paths; indeed, there must be in any 
line segment of length more than 1 . 



Lemma 4. Let p, > log(3/2) be given. There are constants B and Ti (with B as 
in Lemma\^) with the following property. Whenever we start X^(-B) at a point 
of type 'h' in the set Sq = [—B,B] x [-B — 1/2,— B] shown in Figure^ the 
expected number of points of generation Ti of type 'h' in the set Si — [-B, B] x 
[B - 1/2, B] is at least 2. 

Proof. Let A, B and T\ be as in Lemma |31 and let (X' t ) denote the genera- 
tions of the branching random walk under consideration. From the definition 
of Xp, the (non-normalized) probability density function describing the x- or 
y- displacement in a single step is positive on the whole real line. Hence the 
probability that X' 2 contains a point of [— A, A] 2 , which depends on (2:0,2/0)1 is 
bounded below by some constant c = c{B) > 0, independent of (xq, yo). Using 
part (i) of Lemma [3] and induction on k, it follows that 

mx' kTl+2 n[-A,A] 2 \)>2 k c 

for every k. 

Arguing as at the start of the proof, there is some c' > such that for any 
point (x, y) 6 [—A, A] 2 , if we start a copy of 3C^ at a point (x, y) of type 'h' and 
restrict to [— B,B] 2 , the expected number of points of Si that we obtain two 
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generations later is at least d . It follows that 

mx' kTi+4 ns 1 \)>2 k cc'. 

Choosing k large enough that 2 k cc' > 2, and taking k even so that all points of 
generation T2 = kT\ + 4 have type 'h', the result follows. □ 

The same proof shows that Lemma |4] still holds if, instead of looking for 
points of type 'h' in Si, we look for points of type 'v' in the vertical strip 
S[ = [B — 1/2,5] x [-B, B] just inside the right hand edge of [-B, B] 2 . The 
same is true in the following two corollaries. 

Corollary 5. Let fi > log(3/2) and e > be given, and let B and T2 be as 
in Lemma [^J Then every large enough constant N has the following property. 
Whenever we start X M at N points of So of type 'h' and restrict to [— B,B] 2 , 
the probability that generation T2 contains at least N points of Si with distinct 
x-coordinates is at least 1 — e. 

Proof. Let B and T2 be as in Lemma [H and let Zi be the number of T2th gen- 
eration descendants of the iih starting point that lie in Si. Since we start N 
independent (restricted) branching walks, the random variables Zi are indepen- 
dent; Lemma 0] shows that each has expectation at least 2. Although the Z^ are 
not identically distributed, as T 2 is constant we have a common upper bound 
(of order (2e p — 2) 2T2 ) on ^(Zf) for any i. Hence Z = J^iLi Zf has expectation 
at least 2N and variance O(N). Consequently P(Z > N) > 1 — e if N is large 
enough. 

Recalling that T2 is even, in going from generation T2 — 1 to T2 we take 
horizontal steps, so the x-coordinates of the points we reach are distinct with 
probability 1. □ 

The supercriticality of the restricted walk and the survival probability bound 
(|5|) have the following consequence. 

Corollary 6. Let fi > log(3/2) and e > be fixed, let B be as above, and 
let N be a constant. There is a constant T 3 such that, if (B) is started ( as 
usual) with two particles at (0, 0), one of type 'h' and one of type 'v', then with 
probability at least (f>(fx) — 4e either generation T3 or generation T3 + 1 contains 
at least N points of Si with type 'h' with distinct x-coordinates. 

Proof. This follows easily from Lemma [^ii), using arguments similar to those 
in the proofs of Lemma [4] and Corollary We omit the details. Note that 
we need to consider two generations, T3 and T3 + 1, as it may be that only 
one of the initial particles survives, in which case, after a while, either only odd 
generations or only even generations contain particles of type 'h', and we cannot 
say in advance which it will be. □ 

We are now ready to complete the proof of Theorem [1] 
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Proof of Theorem [7J We have already proved the upper bounds on component 
sizes; it remains to prove the lower bound implicit in l[2]). i.e., to show that 

liminf6»(a>,A/w)(A/w)- 1 > 0(A) (6) 

UJ — >oo 

whenever A > log(3/2). As remarked above, it follows by standard results that 
if lo is large enough (so 6(u>, X/lo) > 0), then with probability 1 there is a unique 
infinite component. 

Let A > log(3/2) be fixed, and set \i — A for compatibility with our branching 
process notation above. (Often, at this point in the proof one would fix an 
arbitrary /i < A, but here the usual 'elbow room' turns out not to be needed.) 

Fix £ > with lOe < 4>(n), and let B, N, T2 and T3 be defined as in the 
corollaries above. As before, we condition on (0, 0) being a vertex of G WiP , and 
let (Y t ) denote the neighbourhood exploration process in G WiP , p = h/lo. We 
have already seen that, for any fixed k, the first k generations of the rescaled 
process (Y t /oj) (defined by Y t /oj — {(x/u},y/uj) : {x.y) £ Y t }) may be coupled 
with the branching random walk X 2 — (X t ) so that with probability 1 — o(l) 
as lo — > 00 the first k generations agree up to displacements of individual points 
by up to k/to. 

This coupling result holds also if we explore only the subgraph of G WjP in- 
duced by vertices in [— loB, ujB] 2 , and replace X 2 by its restriction X 2 (£>) to 
[— B, B} 2 . Combining this observation with Corollary [5] above, we see that if lo 
is large enough then, with probability at least 0(/i) — 5e, exploring G W)P from 
(0, 0) within [— loB, loB} 2 , in either T3 or T3 + 1 steps we reach at least N vertices 
«!,..., Kjv of cuSi = [— loB,loB] x [lo(B — 1/2), ojB], with the last step to each 
vertex being horizontal. Note that we may assume that Vx, . . . , vn have distinct 
x-coordinates (since the corresponding particles of do). 

For (a, b) £ Z 2 , let Q a ,b be the square 

Q a fi = ((2a - 1)ujB, (2a + 1)ujB) x ((26 - l)uB, (26 + l)uB) . 

Since we exclude the boundaries, the squares Q a ,b are disjoint. Since N is 
constant (i.e., independent of lo), the coupling above extends to explorations 
started at any TV points of Z 2 , provided that there are no clashes at the first step, 
i.e., that we do not attempt to explore vertically from two points on the same 
vertical line, or horizontally from two points on the same horizontal line. In what 
follows, we can always assume that this proviso is satisfied (from the 'distinct x- 
coordinates' conclusions of Corollaries [5] and [6l and the corresponding 'distinct 
^-coordinates' conclusions of the variants where we end with points near the 
right of [— B, B} 2 ). We shall not comment on this annoying technicality further. 

Suppose that, without testing points in Q a ,b, we have found N vertices 
Vi, . . . , vn of G WjP lying just below Q a ,b, i-e., in the strip of width lo/2 bordering 
Qa,b from below. Then, from Corollary O provided lo is large enough, with 
probability at least 1 — 2s we may find paths in G WjP from {vi} to N points Wi 
of Q a ,b just below the top of Q a ,b, using only vertices in Q a .b- Similarly (using 
the variant of Corollary[5]with S[ in place of Si), with probability at least 1 — 2e 
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we may find such paths to N points just inside the right-hand side of Q a ,b- By 
symmetry, the same conclusions hold if we start with N points Vi just to the 
left of Q a . b . 
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Figure 4: A schematic depiction of the proof that percolation occurs. Roughly 
speaking, whenever there is an arrow joining two shaded regions, from TV vertices 
in the first shaded region we are very likely to find N in the second. We start 
by showing that from vq we find N vertices in the shaded region of Q0.-1 with 
probability close to 

To complete the proof, it is convenient to explore the component Co of G^^ 
containing the point vq = (0, —2Buj) S Qo,-ii hi a manner shown schematically 
in Figure [4] Recalling that vq is present in G WiP with probability p, by Corol- 
lary[6l with unconditional probability at least — 5e) we find paths within 

G u ,p (~l Qo,-i from Vq to N points just below the top of Qo,-i- Let us denote 
this event by E. 

Let us now explore within the squares Q a ,b, a, 6 > 0, working upwards in 
layers, i.e., in increasing order of a + b. When we come to a particular square 
Q a ,b, if we have found N vertices of G w , p near the top of the square Q a ,b-i 
below, or N vertices near the right of the square Q a —i,b> then with probability 
at least 1 — 4e, testing only points inside Q a ,b, we find paths from these vertices 
to N points near the right of Q a .b, and to N points near the top of Q a ,b- If this 
happens, we declare the point (a, b) 6 N 2 to be open. If we did not reach ./V 
vertices of G UiP near the top of the square below or the right of the square to 
the left, we declare (a, b) to be open regardless of what happens inside Q a .b- 

Because each exploration is confined to its own square, each (a, b) is declared 
open with conditional probability at least 1 — 4e. Thus the distribution of open 
points (a, b) G N 2 stochastically dominates a distribution in which each point is 
open independently with probability 1 — 4s. 

Suppose that E holds and that there is an infinite oriented path P in N 2 
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starting at (0,0), in which every vertex is open and every edge goes up or to 
the right. Then vo lies in an infinite component of G W;P . It is easy to check 
that the probability of the existence of such a path P tends to 1 as e — > 0. (For 
example, one can bound the expected number of 'blocking cycles' as in the proof 
of Lemma 6 in [7], and show that this expectation tends to as e — » 0.) Thus 
there is a function /(e) with f(e) — ► as e — > such that 

> (1 - f(e))P(E) > (1 - /( e )(/i/ w )(0(/i) - 5e) 

whenever a; is large enough. Since e > was arbitrary, this proves ©. □ 

We now turn to the finite graphs G n ,u>,p and G' nw p . 

Proof of Theorem^ Let v = {x,y), I < x,y < n, be any potential vertex of 
Gn,u,p, and let C„ and C£ denote the components of u in the graphs G n ,u,p and 
G WiP . We may regard G n ,o),p as a subgraph of Gr WjP , so 

P(|a| >fc) <P(|C;| >fc). 

If A < log(3/2) is constant and lo is large enough, then by Theorem [1] the 
second probability above is at most (A/w)e~ aA< - fc-1 - ) for every k > 1, where 
a\ > depends on A only. Setting A = 3/ax, say, taking = Alogn, and 
recalling that there are n 2 choices for u, it follows that with probability 1 — o(l) 
the graph G n>Ut x/uj contains no components of size larger than j41ogn, whenever 
lo is large enough. 

Turning to G' n u p , although this cannot be seen as a subgraph of G WjP , 
for any given vertex v of G' n u p we may couple the neighbourhood exploration 
processes in G' n u p and in G WjP started at v so that the latter dominates. (Al- 
ternatively, note that the branching process domination argument we gave in 
Guj,p applies just as well in G' nup .) Thus the bound above holds for G' n uj p also. 
Since G niU} , P is a subgraph of G' n u _, the first statement in Theorem [2] follows. 

It remains to show that in the supercritical case, if A > log(3/2) is constant 
and lo — uj(n) with u> — > oo and n/u) — > oo, then G n ^ \/ui and G' n u contain 
a unique giant component of the 'expected' size, namely (4>(X) + o(l))n\/u>. As 
the arguments are mostly rather standard, we shall only sketch them. 

The local coupling arguments above show that we have the expected number 
of vertices in 'small' components: it remains only to show that whp there is a 
component of size at least (0(A) + o{l))n\/uj. Fix e > and set A' = A — e; we 
assume that e is such that A — 2e > log(3/2). It is also convenient to assume 
that A < 1. We start by analyzing components in G n ul y /ui using the extra 
vertices of G„ iWi a/w for later sprinkling. 

We know from Lemma |3] that if G is a large enough constant (depending on 
A' and e), then the branching random walk £\r_ e (C) restricted to [— G, G] 2 is 
supercritical, and indeed survives with probability at least <f>(\' — e) — e whenever 
we start it with two particles of type 'h' and V at a point not too close to the 
boundary of [— C, G] 2 . Let G be the subgraph of G nM .\' induced by vertices 
in [0, 2wG] 2 . When we explore the component of a vertex of G not too close 
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to the boundary, there are two sources of error in coupling this exploration to 
the corresponding (restricted) branching random walk: the first comes from ap- 
proximating the discrete distribution of points on a horizontal or vertical line by 
a Poisson process, which gives rise to an error probability of 0(X' /uj) = 0(1/ uj) 
at each step. The second comes from omitting tests for points on lines already 
claimed by other points; we may account for this by reducing the branching 
process parameter slightly. If we have found m points so far, each claims only 
one line (apart from the first). It follows that for some L = Q(eu), we can 
couple our exploration to dominate 3£\'- e (C) as long as m < L, with an error 
probability of at most A' /u> in each step (coming from the approximation of a 
binomial distribution by a Poisson), and thus a total error probability that is 
0(e). 

Let N be the number of vertices of G that are in components of size at 
least L. From the remarks above it follows that EN > xN , where x — 4>(X' — 
e) — 0(e) — (f)(X) — 0(e), and iV = (2Cuo) 2 X' /ui is the expected number of 
vertices of G. Furthermore, the number N is concentrated in the sense that 
E\N — EiV = O(eNo) (at least for large uj); this can be seen, for example, by 
comparing with Nk, the number of vertices of G in components of order at least 
k, for some fixed but large k, and noting that starting our exploration from 
two different vertices and using a corresponding upper bound, we find that the 
variance of Nk is o(N 2 ). Hence, N > (x — e 1//2 )A^ with probability at least 
1- 0( e V2). 

Sprinkling extra vertices with density (A — X')/u), it is easy to check that, 
whp, all components in G of order at least L join up to form a single component, 
which thus has size at least N. One crude argument is as follows: suppose that 
before sprinkling we have components C\ and C2 each with at least L = Q(u>) 
vertices. Let A denote the set of ^-coordinates of vertices of C\, and B the 
set of y-coordinates of vertices of C2 ■ Since the expected number of vertices of 
G an any horizontal or vertical line in [0, 2ujC} 2 is 0(1), whp no line contains 
more than 0(loguj) vertices, so we may assume that \A\, \B\ = fl(u)/ logui), 
say. Let S be the set of sprinkled vertices in A x B. Then |5| is concentrated 
about its mean, which is of order at least w/(logw) 2 , so we may assume that 
\S\ > uj/(\oguj) 3 , say. Using once more the fact that there aren't too many 
vertices on any line, we may find a subset S' in which all coordinates are distinct, 
with \S'\ — fl(w I (log w) 5 ), say. Finally, for each (x,y) E S', we look for a 
vertical path of sprinkled vertices (outside Ax B) joining (x, y) to the (nearest) 
vertex of C\ with the same x-coordinate, and a horizontal path to Ci- We find 
such paths if each of at most 8C intervals of length u/2 contains a sprinkled 
vertex, an event whose probability is bounded away from zero. These events 
are independent for different (x,y) € S', so with very high probability at least 
one such pair of paths is present. 

Returning to G ntU ,,\/ u , or G' nojX ^, we may cover the vertex set of this 
graph by squares of side- length 2Cu>, overlapping in regions of width Cu>, say. 
The argument above shows that within each square, we almost certainly find 
a component of the right size, i.e., containing at least a proportion x — e 1 / 2 
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of the vertices. The sprinkling above also shows that the giant components 
of overlapping squares are very likely to meet. Considering either mixed or 1- 
independent percolation on a K by K grid in which each vertex/bond is open 
with probability at least p, by estimating the expected number of blocking 
cycles it is easy to check that the probability that any two vertices are joined 
is at least f(p), for some function f(p) tending to 1 as p — ► 1, independent of 
K. Comparison with such a model shows that as n,n/ 'u> — > oo, whp the large 
components in almost all squares are linked up, giving a giant component in 
G n ,u,\/u or G' n u containing at least a proportion </>(A) — 5(e) of all vertices, 
for some function 5(e) — > as e — > 0. □ 

Note that the condition n/co — » oo in the second part of Theorem [2] is 
essential. Indeed, if n/u) is bounded, say n < Ceo, then the approximation 
of the component exploration by the branching random walk X\ breaks down 
even during the first step: this branching random walk has probability 0(1) of 
starting by moving to a point at distance at least C (corresponding to distance 
at least Cuj > n in the graph) from the initial point. This condition was not 
needed in the much weaker result of Frieze, Kleinberg, Ravi and Debany [9], but 
this is not surprising. The very simple analysis in [9] ignored these long-range 
steps of the branching random walk (as well as other complications), which is 
why it does not give the correct limiting size of the giant component, or indeed 
the asymptotic critical probability. 

2.3 The case n = Q(u) 

In this subsection we briefly consider the emergence of the giant component in 
Gn }UJ} p or G' n u p when ui — > oo and n = Q(u>). We assume throughout that 
n ~ Cu> for some constant C > 0. As before, the relevant normalization is to 
take p — p(n) — \/u>, with A constant. 

As noted above, Theorem [5] does not apply in this setting. However, our 
methods still allow us to obtain the critical point, and indeed the size of the 
giant component, in terms of the solutions to certain non-linear equations. As 
this problem is of rather specialized interest, we shall only outline the results 
and arguments. 

Let V\ be a Poisson process on K with intensity A. For d > 1, let g\(d) 
be the probability that within the interval [0, d] there is an interval of length 
1 containing no points of V, and let r\(d) = 1 — g\(d). For < d < 1 we 
set g\(d) — and r\(d) = 1. Given that both and d are points of V\, then 
r\(d) is the probability that they are joined in the graph with vertex set V\ 
and edges between all points at distance at most 1. The function r\(d) will be 
key to our analysis. (This function r\(d) has been studied in other contexts. 
It is the probability that the largest gap between Po(Ad) points uniformly and 
independently distributed in [0,1] is at most 1/d; this largest gap is also a 
version of a scan statistic; see Glaz, Naus and Wallenstein [11] . for example. 
An explicit, but rather complicated, expression for r\(d) is easily derived from 
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the formula given by Stevens [II] for the case of the largest gap between a fixed 
number of random points on a circle.) 

Let us start with the simpler toroidal case. Let Pa,c be a Poisson process 
of intensity A on the circle Tq = K/(CZ) with circumference C. Adding in the 
extra point 0, let Z = Z(X, C) denote (the distribution of) the number of other 
points of V\.c that may be reached from in steps of size at most 1. Let jEa.c 
be the branching process in which the number of children of initial particle has 
the distribution of the sum of two independent copies of Z, and the number 
of children of each later particle has the distribution of Z. Approximating the 
component exploration in G' n w p by the corresponding branching random walk 
on the torus, the proof of Theorem [5] adapts easily to show that if p = X/lu and 
n ~ Cu>, then the asymptotic fraction of vertices of G' n u „ that are in the giant 
component is simply the survival probability p(X, C) of £a,c- ^ n particular, the 
critical value A c = A C (C) of A may be found by solving E(Z(A, C)) = 1. 

If the point x £ Tc is present, then it may be reached from with probability 
1 — g\(x)g\(C — x), so A c is given by the solution to 

c 

(l - g\{x)g\(C -x))dx = l. 

The case of the grid is more interesting, but also more complicated. In this 
case we must consider the restriction £a,c °f our original branching random 
walk X\ to a square [0,C] 2 . Previously, we used a similar restriction with C 
large as a tool in our analysis. Here we are forced to analyze the restriction 
itself, and have no control over C. 

Fixing A and C, for < x, y < C, let ph(x,y), p v {x,y) and p{x,y) denote 
respectively the survival probabilities in the cases where we start with one par- 
ticle of type 'h' at (cc,j/), one particle of type 'v' at (x,y), and one particle of 
each type at (x,y). By symmetry we have ph(x,y) = p v (y,x). Also, 

P(x, y) = 1 - (1 - Ph(x, y))(l - p v (x, y)) = p h (x, y) + p v (x, y) - p h (x, y)p v (x, y). 

Starting with two particles, one of each type, at a point (x, y) chosen uniformly 
from [0,C] 2 , the survival probability p satisfies 

pC pC pC pC pC pC 

P = I / p(x,y)=2 / ph(x,y)- / p h {x,y)p h (y,x), 
Jo Jo Jo Jo Jo Jo 

so to understand p it suffices to understand p^ix^y). 

It is not hard to convince oneself that the proof of Theorem[2]can be modified 
to show that if p = X/w and n ~ Cuj, then the asymptotic fraction of vertices of 
G n . u .p in the largest component is simply p = p(X, C). Filling in the details may 
well require considerable work, however: for example, to make the sprinkling 
argument work one needs to show that p(A, C) is continuous in A. Let us omit 
the details and turn instead to the study of p(X, C). 

It is easy to see that one can express Ph(x, y) as the maximum solution to a 
certain non-linear integral equation; as the details are not very illuminating, we 
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ignore the size of the giant component and consider only the question of when 
it emerges, i.e., when p(X,C) > 0. At this point, it is tempting to claim that 
arguments of the type given in [5], for example, show that p(X,C) > if and 
only if the corresponding linearized operator has norm strictly greater than 1. 
However, there are some complications: in [5], we assumed that the linear oper- 
ator was symmetric and compact; here it is neither! These complications seem 
to be connected to the interchange of x and y coordinates, or, putting it an- 
other way, to the one-dimensional action of the operator on the two-dimensional 
space [0,C] 2 . It turns out that one can simplify things greatly by considering 
two steps of the branching process at once. 

Starting with a point at (x,y) (of type 'h', say, although it doesn't matter), 
and taking two steps in our restricted branching random walk, the expected 
number of points in [x' + dx'\ x [y 1 + dy'} that we reach is easily seen to be 
X 2 r\{\x — x'\)r\(\y — y'\) dx' dy' . Let Ti be the corresponding integral operator 



The operator T 2 is very well behaved: the kernel is clearly symmetric (with 
respect to swapping (x,y) and (x',y')) and, since < r < 1, the operator is 
compact. Thus arguments of the type given [5] do show that p(X, C) > if and 



At this point, we can return to one dimension: the operator T 2 has a unique 
(up to normalization) eigenfunction 4> with eigenvalue the norm of T 2 . Let T be 
the operator on L 2 ([0, C]) defined by 



From the form of T 2 , it is not hard to check that cf>(x,y) = ip(x)ip(y), where 
tp is the eigenfunction of T with maximum eigenvalue. In particular, ||T2|| = 
|T|| 2 . Although we only sketched the details, one has p(X, C) > if and only 
if || T 2 1| > 1, so in this setting the critical value of A is given by the solution to 
||T|| = 1, noting that ||T|| depends both on C and on A. 

Unfortunately it seems unlikely that one can find ||T|| explicitly; this is to be 
expected. The conclusion is that we find the same sort of connection between 
the giant component in a certain random graph and the norm of an integral 
operator on L 2 ([Q, C)) as seen in [5]. This is especially interesting as one cannot 
directly apply the results of [5]. Indeed, here, the approximating branching 
process/random walk is not of the multi-type Poisson form considered there. 

3 Discussion and variants 

The line-of-sight model G UiP is unusual in various ways. There are several 
random graph models in which one can determine the exact threshold for the 



on L 2 ([0,C] 2 ), so 




only if ||T 2 || > 1. 
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emergence of a giant component in terms of the survival probability of a suitable 
a branching process - the simplest example is G(n,p); a very general inhomo- 
geneous model with this property is that of [5]. In lattice percolation models, 
on the other hand, as a rule one cannot determine the exact critical probabil- 
ity. Penrose [T3] defined a 'spread-out' percolation model as follows: let S be a 
fixed bounded symmetric set in IR 2 . Given parameters u> and p, for every pair 
of points v, w of Z 2 , join them by an edge with probability p if v — w G loS, 
independently of all other such pairs. (Penrose's model is somewhat more gen- 
eral, but this is the essence.) For fixed u>, this is a lattice percolation model, 
and the critical probability cannot be found exactly. However, Penrose showed 
that as oj — ► oo, the asymptotic form of the critical probability can be found - 
the 'critical expected degree' tends to 1. In [6], it was shown that this result 
follows easily from those of [5]. Given the homogeneous nature of the model, 
the critical expected degree tending to 1 shows that asymptotically, cycles do 
not matter; this is far from the case in the line-of-sight model. 

In terms of the critical probability, the line-of-sight model G^p behaves dif- 
ferently from Penrose's model: even when to is large, G WjP contains many short 
cycles, so the critical expected degree, 41og(3/2), is larger than 1. There are 
two important differences between the models. Firstly, G^.p is a site percola- 
tion model (vertices are selected at random), and Penrose's model is a bond 
percolation one (edges are selected at random). Secondly, in G WiP the geometric 
condition for joining vertices involves scaling a set S C I 2 of measure zero. 
We shall examine the effect of these differences separately, by considering two 
variants of G^.p. 

There is a natural bond percolation variant of G WiP . Given u> and p, let G^p 
be the graph on Z 2 defined as follows. For every pair of vertices v, w such that 
the line segment vw has length at most to and is horizontal or vertical, join v 
and w with probability p, independently of all other such pairs. We consider 
the limit w — > oo with p = X/lo, where A is constant. In this model the density 
of short cycles is very low. Thus, as we explore the component of a vertex in the 
usual way, the number of vertices reached at each step may be approximated 
by a Galton-Watson branching process in which each particle has a Poisson 
number of children with mean 4A. Moreover, the location of these vertices 
may be approximated by a corresponding branching random walk Xi\. The 
arguments above for G^ iP carry over to this setting, showing that the critical 
probability p c {oj) is asymptotically l/(4w), and that if p = X/ui and cu — > oo then 
the infinite component contains a fraction ^(4A) + o(l) of the vertices, where 
4>(/J.) is the survival probability of the branching process underlying X M , i.e., the 
maximal solution to 0(/i) = 1 — e _A "^^. One can also obtain corresponding 
results for finite graphs, showing that when uj — > oo, they behave much like the 
Erdos-Renyi (or Gilbert) model G(n,4A/n). 

Turning to our second variant, let C £ be the cross with (four) arms of length 
1 and width 2e, so 

C s = [-1,1] X [-£,£] U [-£,£] X [-1,1]. 
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Let Gt^p.E be the random graph defined as follows: define a graph G on I? by 
joining two vertices v and w if v — w £ luC £ . Then select vertices independently 
with probability p to form G^.p.e- When e — 0, this graph is exactly G u>p . 
For e > 0, the area of wC e is proportional to to 2 , so we consider the limit with 
p = X/lu 2 and ui — > oo. In this context it is natural to rescale the vertex set, 
selecting vertices of (l/w)Z 2 with probability X/lu 2 and joining them if their 
difference lies in C e . In the limit, the selected vertices form a Poisson process 
on R 2 with intensity A. This gives us another natural random graph model, 
G\ :S . The vertex set of G\ tE is a Poisson process V\ on R 2 with intensity A, and 
two vertices v and w are joined if and only if v — w G C £ . It is not hard to see 
that as to — > oo with e fixed, to 2 times the critical probability for percolation in 
Guj,p,e approaches the critical density A c for percolation in G\_ e . 

The random graph Ga, e we have just defined is a special case of a percolation 
model introduced by Gilbert [lOj in 1961, Gilbert's disc model. To define this, 
let S be a symmetric (in the sense S = —5) set in R 2 with Lebesgue measure 
< area(S') < oo, and let G(A, S) be the random graph whose vertex set is the 
Poisson process V\, in which vertices are joined if their difference lies in S, so 
G\. £ — G(A, C e ). Gilbert asked the following question: given S, what is the 
critical value A c = A c (5) above which G(A C , S) contains an infinite component? 
Since G(X,S) and G(X/c 2 ,cS) have the same distribution as graphs, one often 
works with the 'critical expected degree', or critical area a c (S) — X c (S)&ie&(S) , 
instead. (The term critical area is used because of the standard normalization 
A=l.) ^ 

As with many such percolation questions, it seems impossible to determine 
a c (S) exactly. The most studied cases are when S is a disc D, or when S is 
a square [— 1, l] 2 . The best rigorous bounds for a c (D), the bounds 2.184 < 
Qc(-D) < 10.588 proved by Hall [T2], are not much better than Gilbert's original 
bounds; although these bounds are far apart, it seems hard to improve them 
significantly. In practice, however, a c (D) is known quite precisely; Balister, 
Bollobas and Walters [3] proved that [4.508, 4.515] is a 99.99% confidence interval 
for a c (D), and [4.392, 4.398] for a c ([-l, l] 2 ). See [7, Section 8.1] for more details, 
as well as many references to heuristic bounds. 

Recall that we introduced G\, e — G(X,C £ ) as a variant of G w , p . In this 
context it is natural to take the limit e — > 0, and ask whether the critical 
average degree converges to the corresponding value for G UjP . Since the area of 
the cross C £ is 8e — 4e 2 , one might expect that 

8eA c (G £ )~41og(3/2) (7) 

as e — > 0. We shall see that this is not the case. In fact, rather surprisingly, we 
can describe the limit of eA c (G e ) in terms of the Gilbert model with S a square! 

Given A > 0, let G be the graph G(A, [— 1,1] 2 ) conditioned on the origin 
being a vertex. In other words, G is the graph on V\ U {(0, 0)} in which we join 
two vertices if their distance is at most 1. Let 

/ 2 (A)=E(|Go|-l), 
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where Co is the component of G containing (0,0), and [Co is the number of 
vertices in Co- Note that we subtract 1 to avoid counting the 'extra' vertex at 
the origin. Later we shall consider an analogous quantity defined in terms of a 
d-dimensional form of Gilbert's model; this is the reason for the notation f 2 . 

Theorem 7. The critical densities A c (C e ) satisfy e\ c (C E ) — > fi as e — > 7 where 
\x is the unique solution to f 2 (p-) = 1- 

Proof. Note first that the definition of \x makes sense: we have /2(C)) = 0, while 
/2(A) diverges at A c ([— 1, l] 2 ) < 00. Furthermore, up to this point /2(A) is 
increasing and continuous; hence there is a unique solution to fiiji) = 1. 

Fix A > and e > 0; we shall take £ _1 A as our density parameter. 

Let us condition on u = (0,0) being a vertex of Gw £i£ = G(A/e, C e ), and 
explore the component of u in this graph. As in our study of G WjP , we break 
the exploration down into steps; in each step, we explore either horizontally or 
vertically from a vertex v. As before, we set Yq — {vq}, explore horizontally and 
vertically from vo, writing Y\ for the set of vertices reached, and then explore 
horizontally from those vertices of Y t reached during vertical explorations, and 
vice versa. We write Y t +\ for the new vertices found during such explorations, 
so the vertex set of the component containing i>o is the disjoint union of the sets 
Y t , t > 0. This time, however, the horizontal and vertical explorations are a 
little more complicated. 

Let us say that an edge vw of G\/ e , e is horizontal if v — w £ [— 1, 1] x [— e, e]. 
The definition of vertical edges is analogous. Note that some edges are both 
horizontal and vertical, but the proportion of such edges tends to as e — > 0. 

To explore horizontally from v, let H v consist of all new (not reached in 
previous explorations) vertices v' of G\/ Sie joined to v by horizontal edges, 
together with all new vertices v" joined to such vertices v' by horizontal edges, 
and so on. In other words, writing Gj, ££ for the subgraph of G\/ e<e formed 
by all horizontal edges, H v U {v} is simply the component containing v in the 
subgraph of G\j ee induced by the v and the new vertices. In particular, if 
this is our first exploration, then v = v = (0,0), and H v U {v} is simply the 
component Co of the origin in Gj, g . 

Now Gj, s is exactly the Gilbert graph G(A/e, [— 1, 1] x [—e, e]), conditioned 
on (0, 0) being a vertex. Scaling vertically by a factor e, the distribution of this 
random graph (as an abstract graph) is identical to that of G(A, [— 1, l] 2 ). Let 
T>2(X) denote the distribution of \Cq\ — 1, where Co is the component of the 
origin in this graph, so /2(A) = E(2?2(A)) by definition. Then \H Vo \ has the 
distribution V 2 (X). 

In later steps of the exploration, the restriction to new vertices ensures that 
the distribution of H v , given the exploration so far, is stochastically dominated 
by that of H Vo , i.e., by T>2(X). The same holds for vertical explorations. It fol- 
lows that the sequence \ Y \, |Yi |, . . . is stochastically dominated by \Xo\, \X\\, . . ., 
where (X t ) is a branching process in which the number of children of each par- 
ticle except the initial one has the distribution T>2(\), and different particles 
have children independently. The number of children of the initial particle is 
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distributed as the sum of two independent copies of T>2(\), since we explore 
horizontally and vertically from Vq. 

Since E(X>2(A)) = /2(A), this shows that G\/ e e has no giant component 
when X < fi, where fi satisfies f2(t L ) = 1- Hence, e\ c (C e ) < fx. 

It remains to show that if A > /1 is constant, then for e small enough we have 
percolation in the graph G\/ ££ . The argument is very similar to that for G WjP , so 
we give only the briefest outline, emphasizing the differences. In doing so we may 
assume that A < A c ([— 1, l] 2 ): if A > A c ([— 1, l] 2 ) then the horizontal subgraph 
of G\/ e e has the same distribution as a graph as the supercritical Gilbert graph 
G(A, [—1, l] 2 ), so this subgraph already contains an infinite component. 

When exploring horizontally in G\i ejS , there is a priori no bound on how 
far we may drift vertically. To deal with this, let us choose a large constant A 
(depending only on A), and limit each horizontal or vertical exploration to A 
steps. Let X>2 W be the distribution of the number of points reached from 
uo in this restricted horizontal exploration. Since the horizontal subgraph of 
G\/e,e is subcritical, the distribution T> 2 A \X) converges to X>2(A) as A — > 00. In 

particular, since A > /1, we may choose A so that E(2?2^(A)) > !• 

With this modified exploration, the coupling argument used in G UjP carries 
over to show that we may couple the first O(l) steps of the exploration with 
the appropriate branching process with an error probability of 0(e). The main 
difference is that rather than claiming lines, new vertices we find claim hori- 
zontal or vertical strips of width 2Ae. As before, we can in fact couple with 
a supercritical branching random walk, and use oriented percolation to find an 
infinite component. □ 

Let A e be the annulus centered on the origin with outer radius 1 and inner 
radius 1 — e. Independently, Franceschetti, Booth, Cook, Meester and Bruck [8], 
and Balister, Bollobas and Walters [3J showed that lim e ^oac(^£) = 1- This 
result is similar in spirit to the results here: one shows local approximation by a 
branching process (in this case, a very simple process, since the graph contains 
few short cycles) , and then must work to deduce percolation in the supercritical 
case. 

Let S e be the 'square annulus' [—1, l] 2 \ [— (1 — e), 1 — e] 2 . One might expect 
lim e ^o a c (S e ) to also equal 1, but this is not true: in [3J the bound a c (S e ) > 
1.014 is proved for every e > 0. Our methods here show that a c (S e ) docs 
converge as e — > 0, and give a description of the limiting value. 

Let F = [—1,1] 2 x {-l,l},soFcl 2 xZ consists of two opposite faces of 
the cube [—1, l] 3 . We define a variant G(A, F) of Gilbert's model as follows: the 
vertex set is a Poisson process of intensity A on R 2 x Z, with the origin added, 
and two vertices v and w are joined if and only if v — w € F. Let Go be the 
component of the origin in G(A, F), and let /(A) = E(|Go| — 1). 

Theorem 8. The critical area a c (S £ ) tends to 2/j as e — > 0, where /i is the 
unique solution to /(/i/8) = 1. 
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Note that the expected degree in G(A, F) is 8A; thus the limiting expected 
area or degree in Theorem[8]is twice the expected degree at which E(|Cq| — 1) = 1 
holds in G(A, F). In this sense, Theorem [5] is analogous to Theorem [7] 

Proof. The proof is very similar to that of Theorem [71 so we omit the details. 
The key observation is as follows: call an edge vw of G(A, S e ) 'horizontal' if v — w 
lies in one of the vertical sides of the annulus, i.e., if ±(v — w) G [1 — e, 1] x [—1, 1], 
and define vertical edges similarly. As e — > 0, it is easy to check that almost 
all short cycles in G(X/e,S £ ) consist entirely of horizontal edges, or entirely 
of vertical edges. As before, we break down the neighbourhood exploration 
process in G(A/e, S £ ) into horizontal and vertical explorations. In a horizontal 
exploration, say, we look for all new vertices that can be reached by horizontal 
edges. As before, the asymptotic condition for criticality is that the expected 
number of new vertices found in a single horizontal (or vertical) exploration is 
1. Rescaling horizontally by 2/s, a horizontal exploration may be coupled with 
an exploration of the component of the origin in the Gilbert model G(A/2, T e ), 
where T e consists of two squares of side 2 centered at (±(2e _1 — 1), 0). If we 
take k < 1/e steps from the origin, with the displacement of each step in T e , 
then we can tell from our final position how many steps were to the right and 
how many to the left. It follows that for the first 1/e steps, the explorations 
in G(X/2,T E ) and in G(A/2,F) may be regarded as identical. It follows that 
eX c (S E ) ~ 2(/x/8), where /(/x/8) = 1. Since the area of S E is asymptotically 8e, 
the result follows. □ 

Since it is much easier to estimate (by simulation) the point at which an 
expectation crosses 1 than the point at which it diverges, Theorem [3] makes it 
much easier to estimate lim e _»o a c (S £ ); a simple simulation suggests that the 
limit is 1.11406 ± 0.00001. 

The models above make just as good sense in any dimension. For d > 1 and 
A > 0, let G = G(A, [—1, l] d ) be the graph whose vertex set consists of a Poisson 
process in IR d of intensity A together with the origin, in which two vertices are 
joined if they are within distance 1. Let 

/ d (A) = E(|G |)-l, 

where Go is the component of G containing the origin, so fd(X) = oo if A > A c , 
where A c is the critical density for G(A, [—1, l] d ). In general, one cannot hope 
to evaluate fd(X) exactly. However, when d = 1, Go — 1 has exactly the 
distribution defined earlier, so /i(A) = 2(e A — 1), and /i(log(3/2)) = 1. 
It is easy to check that / 2 (A/2) > /i(A) for any A > 0. (This is the natural 
comparison, as the expected degree of a vertex of the relevant graph is 2 d X.) It 
follows that the quantity fi defined in TheoremUJis strictly less than log(3/2)/2, 
so Theorem [7J shows that our 'guess' © does not hold. (A quick simulation 
suggests that [i is around 0.177635.) In other words, the critical expected degree 
in models defined using a cross with very thin arms does not approach that of a 
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model defined using a cross whose arms have width 0. (Compare 8 x 0.177635 = 
1.421 . . . with 4 log(3/2) = 1.621 . . . .) 

We finish by turning to one final class of variants of G UiP , namely the nat- 
ural generalizations to higher dimensions. Of course, there is more than one 
natural equivalent of a 2-dimensional cross in higher dimensions: a union of line 
segments, or a union of (d — l)-dimensional hypercubes. Fortunately one can 
cover both generalizations in a single definition. 

Given d>2, 1 < r < d — 1, u> >1 and < p < 1, let Gd, r ,ui,p be the random 
graph defined as follows: for the vertex set, select points of Z d independently 
with probability p. Join two vertices v and w if and only if v and w differ in 
at most r coordinates, and differ in each of these coordinates by at most w. If 
d — 2 and r = 1 then we obtain G^p. 

Theorem 9. Let p c {d,r,u>) denote the critical probability p for percolation in 
Gd, r ,uj,p- Then 

Lu r p c (d,r,uj) -» X d>r 

when u> — > oo with d > 2 and 1 < r < d — 1 fixed, where Xd, r * s the unique 
solution to 

((^)-l)/,(A rf , r ) = l. 
Recalling that A (A) = 2(e A - 1), we see that 

so TheoremlHlgenerahzes the asymptotic critical probability result for G WjP given 
by Theorem [T] Note that the function /2OU), which appeared in Theorem [7] in 
the analysis of a certain 2-dimensional graph, appears here in the analysis of 
the very different graphs Gd,2,w,p ; d > 3. 

The proof of Theorem [9] is very similar to that of Theorem [1] so we omit it. 
The main difference is that instead of exploring horizontally or vertically, i.e., 
always exploring within a line, we always explore within an affine subspace of 
dimension r generated by r of the coordinate axes. The factor (f) - 1 above 
appears because there are (^) such subspaces through any given point, and we 
reach any vertex other than the initial one by exploring within one of them. 
This time, when we find a vertex v of Gd. r .ui,p by exploring within a certain 
subspace, this vertex immediately claims the other ( ) — 1 subspaces in which 
it lies. As before, we omit testing points w in two (or more) claimed subspaces. 
When exploring from v in a subspace S, any such point w lies in a subspace S' 
through a vertex v' of Co that we have already found, with S' and S unequal and 
therefore not parallel. Since S and S' intersect in a space of lower dimension, 
in the early stages of our exploration almost all points of S are claimed only by 
v, and the coupling goes through as before. 
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Just as for G WiP , we could 'thicken' the generalized crosses defining Gd, r ,uj,p 
by e and let e — > 0, obtaining a generalization of Theorem [7] In this case the 

limit jj, in Theorem [7] is replaced by the solution to f(^) — l^j /d(/i) = 1. 
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